Vorticity affects the stability of neutron stars 
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The spin rate of neutron stars at a given temperature T is constrained by the interplay between 
gravitational-radiation instabilities and viscous damping. Navier-Stokes theory has been used to 
calculate the viscous damping timescales and produce a stability curve for r-modes in the (fi, T) 
plane. In Navier-Stokes theory, viscosity is independent of vorticity, but kinetic theory predicts 
a coupling of vorticity to the shear viscosity. We calculate this coupling and show that it can in 
principle significantly modify the stability diagram at lower temperatures. As a result, colder stars 
can remain stable at higher spin rates. 
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It has recently been shown [|h\|2j] that in rotating stars 
the r-modes are unstable via a coupling to gravitational 
radiation. This instability can account for the spin-down 
of hot newly-born neutron stars. Viscosity damps out 
the oscillations, and tends to stabilize them in general. 
The timescales associated with this instability have been 
calculated for uniform density and Newtonian polytropic 
stellar models up to the lowest order of stellar spin rate 

r-mode instability could be crucial in understanding 
the observed spin rates of neutron stars at different tem- 
peratures, and it is also a potentially observable source 
of gravitational radiation It is therefore important 
to improve our understanding of the various factors that 
go into r-mode stability analysis. As examples, the bulk 
viscous timescale has recently been calculated to second 
order in Q and relativistic effects on r-modes have 

recently been considered [|| . Further work remains to be 
done, both in refining the thermo-hydrodynamical model, 
and in investigating relativistic r-modes. In this Letter, 
we consider the possible effects of vorticity on shear vis- 
cosity. These effects are predicted by kinetic theory when 
the unperturbed equilibrium state is rotating. 

In standard Navier-Stokes theory, the angular velocity 
of the fluid has no effect on viscous stress or heat flux, 
which obey the equations 



II = -C© , q a = -tV„T , TT ab = -2rj(j ab , 



(1) 



where IT is the bulk viscous stress and C, is the bulk viscos- 
ity; q a is the heat flux and k is the thermal conductivity; 
ir a b is the shear viscous stress and r\ is the shear viscosity; 
= V ' a v a is the volume expansion rate of the fluid, T 
is the temperature, and a a b = V \ a v b ) is the rate of shear 
(where the angled brackets denote the symmetric trace- 
free part). The fluid vorticity uj a = ^curlu a does not 
enter these equations, even when the equilibrium state is 
rotating. 



On physical grounds, one might expect that rotational 
accelerations can couple with gradients of momentum 
and temperature, so that there could in principle be cou- 
plings of uj a to q a and ir ao . In the case of heat flux, 
qualitative particle dynamics indicates || (p. 34) that 
this coupling does exist as a result of a Coriolis effect, 
which is in some sense analogous to the Hall effect in a 
conductor subject to a magnetic field. The Coriolis effect 
on heat flux is confirmed by molecular dynamics simula- 
tions @. Miiller @ and Israel & Stewart @ showed 
that the Boltzmann equation predicts in general a cou- 
pling of vorticity to heat flux and shear viscous stress. 
The microscopic and self-consistent kinetic approach is in 
contrast to the continuum view, where a phenomenolog- 
ical principle of "frame indifference" is invoked to argue 
against any vorticity coupling. (See |) 10 1^] for further 
discussion.) 

Using the Grad moment method to approximate the 
hydrodynamic regime via kinetic theory, the relations in 
Eq. (0) are modified to f| (Eq. (7.1)) 

n = -c [e + p ti] , (2) 

q a = -K [V a T + T{3\ {q a - UJ ab q b }] , (3) 
TT a b = -2r/ [(Jab + P2 (ab) ~ 2^ C ( Q 7T fc ) c }] , (4) 

where Pa can be evaluated in terms of collision integrals 
for specific gases, an overdot denotes the comoving (La- 
grangian) derivative, and the vorticity tensor is given by 

Uab = V[ Q W b ] = e aoc uj c , 

where square brackets on indices indicate the skew part. 
Navier-Stokes theory is recovered from the Miiller-Israel- 
Stewart theory when Pa = 0. However, kinetic theory 
gives Pa values for simple gases which are definitely not 
zero. Furthermore, if Pa — 0, the equilibrium states 
are unstable and dissipative signals can propagate at un- 
bounded speed OJl3|. 
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The /^-corrections will be very small except if there 
are either high frequency oscillations (pumping up the 
time-derivative terms) or rapid rotation (pumping up the 
vorticity-coupling terms). In the context of rapidly ro- 
tating neutron stars, we expect the vorticity-dissipative 
couplings to dominate the time-derivative terms; this ex- 
pectation is borne out by calculations (see below). The 
vorticity-dissipative couplings will be negligible if the un- 
perturbed equilibrium state is irrotational, i.e., if uj a = 
in the background, so that the coupling terms become 
second-order. However, for fast rotation, uu a ^ in the 
background and the coupling terms make a first-order 
contribution to dissipation. In the words of Israel & 
Stewart p|i : 

"these results will ultimately be of practical interest in 
astrophysical and cosmological situations involving fast 
rotation, strong gravitational fields or rapid fluctuations 
(neutron stars, black hole accretion, early universe), al- 
though it will probably be some time before the state of 
the art in these fields makes such refinements necessary." 
We believe that recent and ongoing developments in ro- 
tating neutron star physics have reached the stage where 
the Miiller-Israel-Stewart theoretical corrections to the 
Navier-Stokes equations need to be examined, and our 
results indicate that the corrections could be important. 

We follow the standard assumption |Q that the heat 
flux may be neglected relative to viscous stresses in cal- 
culating damping timcscales. (Elsewhere we will dis- 
cuss how the vorticity-heat flux coupling could affect 
this standard assumption.) Then the vorticity correc- 
tion to Navier-Stokes theory reduces to the coupling term 
cjC (a 7r 6)c- This term means that the angular momentum 
of the star changes the shear viscosity timescale, and 
we find (for axial r-modes) a correction proportional to 
T~ r fl 2 , where r = 9 for a nonrelativistic fluid and r = 12 
for an ultrarelativistic fluid. 

The evolution of dissipation energy contained in small 
fluctuations is given by 



dE 
~dt 



c 



d"7r ab <57r* 
2r? 



d A x 




(5) 



where (dE/dt) G is the energy flux in gravitational radia- 
tion, oTI = n — n and 5it ab = n ab — n a b) with an overbar 
denoting background quantities. In this case, ft = = 
n ab . The normal modes of the star are damped by dissi- 
pation, and the damping rate can be determined by Eq. 
(||). For a normal mode with time dependence e lro *, the 
energy has time dependence exp[— 2Im(tn)i]. Then by 
Eq. (|^), the characteristic damping time r = 1/Im(ro) 
of the fluid perturbation is given by 



1 dE 
2E dt 



1 



(6) 



where r B , r s , and r G are the bulk viscous, shear viscous, 
and gravitational radiation timescales respectively. 



To evaluate the vorticity-corrected shear viscous 
timescale, we use Eq. (§) in E qs. (|) and @. To lowest 
order 



STTab 



-2r? 



5a a b - 2i?77?7/?2oV a6 + Ar)f3 2 Sa? a uj 



b)< 



where uj a is the background vorticity (the background 
shear vanishes). Then 

5ir ab 5ir* ab = V {Sa ab Sa* ab + 4 7 2 [w 2 5a ab 5a* ab 

+ 4 (Sa ab Sa* ab Lu c Lu c - 6a ca 6a* da uj c u d )] } , 

where 7 = 77/32 - The first term is the usual term in 
Navier-Stokes theory, while the following terms are the 
Miiller-Israel-Stewart corrections. The zu 2 term arises 
from 7T a (, in Eq. (|j), and is negligible relative to the uj 2 
terms which arise from the lu c ( a iTfy c term in Eq. (Q). The 
energy dissipation rate through shear viscosity will be 

--ijn {5a ab 5a* ab - 4 7 2 [^Sa ab Sa* ab 

+ 4 (5a ab 5a* ab oj c uj c - 5a ca 5a* da uo c u d )} } d 3 x. (7) 

In order to proceed further, we need expressions for 
the shear viscosity r\ and the coupling coefficient /?2 ■ For 
the various interactions, r)(p,T) is calculated in p^| , ^6[ , 
where it is shown that electron-electron scattering is more 
important for shear viscosity than other interactions. 
The expression for 77 is given in Jl4| , in good agreement 
with (T5]|l|], as 




9/4 



1.10 x 10 16 



10 14 g/cin 



10 9 K 
T 



;/cms. (8) 



For a Maxwell-Boltzmann gas, the coefficient /?2 is found 
in [^2), but we require the expression for a degenerate 
Fermi gas. This has been found by Olson & Hiscock JTt| ] 
in the case of strong degeneracy: 



/3 2 = 



15n 2 h 3 + 
m^gc 5 (y 2 + 2v) 5 / 2 



O 



kT 



(9) 



where m is the particle mass, g is the spin weight, and 
mc 2 v/kT 3> 1. The dimensionless thermodynamic po- 
tential v = (j) + p)/nm — mc 2 s/kT — 1, where s is the 
specific entropy, is equal to the nonrelativistic chemical 
potential per particle divided by the particle rest energy. 
For a strongly degenerate gas, the nonrelativistic chemi- 
cal potential is proportional to T, so that 

kT 



where a 1 is a dimensionless constant measuring the 
degree of degeneracy. The nonrelativistic regime is ob- 
tained for while the ultrarelativistic case corre- 
sponds to v S> 1. 
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For temperatures below 10 10 K, neutrons in the neu- 
tron star are nonrelativistic, while electrons are ultrarel- 
ativistic [fill. The nonrelativistic limit of /3 2 is 



cm s 2 /g, 



(/3 2 ) NR « 3.16 x 10- 5 (aT)- 5 / 2 
and its ultrarelativistic limit is 

(/3 2 ) UR » 6.45 x 10 15 (aT)- 4 cm s 2 /g 
Using Eqs. (|), © and (|ll|), we have 



1.10 x 1CT 11 



(10) 
(11) 

(12) 
(13) 



In the calculation, we used the same relation for rj in both 
cases, because in the high-density regime (p > 10 14 g/ 
cm 3 ) for both electron-electron scattering and electron- 
neutron scattering, r\ is proportional to T~ 2 , with nearly 
equal proportionality factor pl|. For typical values of 
the temperature, T = 10 9 K, and density, p — 3 x 10 14 




g/cm , we find that j v 



x 10 s, while 7 NR 



a -5/2 x 10 -io s ^ 

We assume that the background is a uniformly rotating 
star, so that the equilibrium fluid velocity is v a = fl(fi a , 
where <~p a is the rotational Killing vector field The 
vorticity vector of the equilibrium state is 



0=£.[cot0,-l,O] . 
2r 



(14) 



The r-modes of rotating barotropic Newtonian stars have 
Eulerian velocity perturbations given by H 



Sv = CRfl Y£ exp(irot) , 



(15) 



where C is an arbitrary constant, R is the unperturbed 
stellar radius, and w — 2m£l/£(£+l). The magnetic-type 
vector spherical harmonics arc defined by 



Y B - 



= V X 



(16) 



The shear of the perturbed star is given by 

S<Tab = V <Q (5u b ) . (17) 

Substituting Eqs. into Eq. (0), we find the 

shear viscosity timescale for I = m: 



1 



(^-1)(2^ + 1) / rfr" dr + Q 2 S e 



(18) 



where QJ 1 = J„ pr u+2 dr. The first term in brackets is 
in agreement with the expression calculated in Q, and 
St is the correction term: 



S f re 16 



+ 



(I -1) (21+1) 
(^ + 1) 2 U ° 

l(l-2)![(2l-l)!!] 2 r(|) 



(£+l)(2£-l)(2£)! T{£-\) 
x [(2£ 3 - 8£ 2 -U- 6)U 2 + 12{£ 3 - £ 2 - £ 
+ 2(A£ 4 -£ 3 - 9£ 2 + 5£+ 1)U 4 ] , 

where U k {T) = R k j 2 w 2l ~ k dr. 

For the £ = 2 modes, Eqs. @ and @ give 

— = 5Q 2 / dr 

Ts Jo 

+ iQitf [80C/ + 93C/ 2 + 54C/ 3 - 42t/ 4 ] . 



1)03 

(19) 



(20) 



For comparison with previous calculations based on 
Navier-Stokes viscosity (see, e.g., (?]]), we use an n = 1 
polytrope with mass M = 1.4M@ and radius R — 12.57 
km to evaluate the integrals in Eq. (p0|). The bulk vis- 
cous and gravitational radiation timescales are unaffected 
by the vorticity correction, and we obtain 



1 



r(0,T) 
1 

+ — 

r s 



1 

f a 
10 9 K 
T 



l + qa 4 



f B ^lO^K J ^ 
'lO 9 !^ 



T 



(21) 



where fl K = ^/irGp, which is | times the Keplerian (mass- 
shedding) frequency, and the vorticity correction factors 
are 
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1.36 x 10~ 23 , 
5.67 x 10- 10 , 



9 nonrcl , 
12 ultrarel. 



(22) 



The standard result (see, e.g., 0) is regained for q = 0, 
with 

f G = -3.26 s , f B = 2.01 x 10 11 s , f s = 2.52 x 10 8 s . 

We note that the contribution from the ~k a b term in Eq. 
(0) to the g-correction is less than 1% of the contribution 
from the uj c ( 7T(,\ C term. 

Now we are able to determine from Eq. ( pl| ) the criti- 
cal angular velocity Cl c , defined by l/r(O c , T) = 0, which 
governs stability of the star: if Q > f2 c , then dissipative 
damping cannot overcome the gravitational radiation- 
driven instability. In Fig. 1 we plot £l c /fl K against tem- 
perature T, showing how the vorticity-viscosity coupling 
affects the standard result (see, e.g., [pj). Electrons are 
assumed to dominate the shear viscosity, and they are 
ultrarelativistic over the range of temperatures. 

It is clear from Fig. 1 that the vorticity correction is 
only appreciable at temperatures T < 10 s K, but that for 
these lower temperatures, the correction can be large, 
especially for smaller a. As the degree of degeneracy 
increases (i.e., with increasing a), the correction is con- 
fined to lower and lower temperatures. The effect of the 
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vorticity-viscosity coupling is to increase the stable re- 
gion, so that cooler stars can spin at higher rates and 
remain stable. This may modify recent results (T^| which 
suggest that r-mode instability could stall the spin-up 
of accreting neutron stars with T > 2 x 10 5 K; if the 
vorticity correction operates, then the stability region is 
increased, so that spin-up could be more effective, espe- 
cially for lower degeneracy parameter a. 

Of course, our analysis is limited by the fact that we 
have followed the standard assumption in viscous stabil- 
ity analysis and ignored superfluid effects that will be- 
come important at lower temperatures (see, e.g., Q). 
Superfluid "friction" effects are thought to prevent /- 
mode instability, and these effects are likely to be rele- 
vant also for r-modes. These effects may strongly alter 
the vorticity correction effect, a subject which is currently 
under investigation. In addition, we have used for our r- 
mode calculations solutions that assume slow rotation. 
Thus the fi/Q K > 0.3 part of Fig. 1 is an extrapolation 
to high spin rates, in common with previous stability dia- 
grams. Recent calculations of r-modes for rapid rotation 
p(| should be used in future calculations of the vortic- 
ity correction. Since /-modes are unstable at high spin 
rate, the effect of the vorticity correction on these modes 
would also be interesting to calculate. 

In conclusion, we have shown that the coupling be- 
tween vorticity and shear viscous stress predicted by ki- 
netic theory can in principle have a significant effect on 
r-mode instability in neutron stars. The Muller-Israel- 
Stewart correction of Navier-Stokes theory predicts that 
colder stars can remain stable at higher spin rates, so 
that accreting spin-up could be protected from r-mode 
instability. 

Finally, we remark that the vorticity correction to heat 
flux, as well as the couplings between heat flux and vis- 
cous stress predicted by kinetic theory flllTJl2"i [but not 
shown in Eqs. could lead to interesting modifi- 

cations of the standard stability curve. In particular, the 
coupling of bulk viscous stress to heat flux could have an 
effect at high temperatures ^ij . 

We are indebted to Nils Andersson and Sharon Morsink 
for crucial comments, and we thank also Marco Bruni, 
Luis Herrera, Mehdi Jahan-Miri, and Josep Triginer for 
helpful discussions. VR was supported by a Royal Society 
grant. 
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FIG. 1. Critical angular velocity versus temperature (n — 1 
polytrope with mass 1.4M© and radius 12.57 km). The stabil- 
ity region is below the curves. The solid curve shows the stan- 
dard result, with no coupling of viscosity to vorticity. Broken 
curves (labelled by the degeneracy parameter a) show how the 
instability region is reduced by the kinetic-theory coupling of 
shear viscosity to vorticity, for an ultra-relativistic degenerate 
Fermi fluid (electron-electron viscosity). 
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